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Transverse-momentum resummation for top-quark pairs at hadron colliders 
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We develop a framework for a systematic resummation of the transverse momentum distribution 
of top-quark pairs produced at hadron colliders based on effective field theory. Compared to Drell- 
Yan and Higgs production, a novel soft function matrix is required to account for the soft gluon 
emissions from the final states. We calculate this soft function at the next-to-leading order, and 
perform the resummation at the next-to-next-to-leading logarithmic accuracy. We compare our 
results with parton shower programs and with the experimental data at the Tevatron and the LHC. 
We also discuss the implications for the top quark charge asymmetry. 
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The top quark is of special importance in the Standard 
Model (SM). Due to its large mass, it couples strongly to 
the Higgs boson, and is crucial to the hierarchy problem. 
New physics (NP) models aiming at solving the hierarchy 
problem often predict top partners which exhibit similar 
properties as the top quark and may decay into it. Possi- 
ble new heavy resonances usually prefer to decay into top 
quark pairs. Therefore, studying the top quarks can on 
one hand help understanding the nature of electroweak 
symmetry breaking, and on the other hand probe NP 
beyond the SM. 

If a heavy resonance decays into a top quark pair, the 
kinematics of the tt system will then carry information 
of the resonance. It is therefore worthwhile to study the 
tt pair as a whole instead of individual top quarks. One 
important example is the invariant mass of the tt pair, 
which is very sensitive to new physics contributions. Pre- 
cision predictions for this distribution has been achieved 
in Besides the invariant mass, another important 
variable is the transverse momentum qx of the tt system, 
which has been recently measured by both the CMS and 
the ATLAS collaborations at the LHC [H, H| ■ One reason 
to study this distribution is that the top quark charge 
asymmetry exhibits intriguing dependence on qx 0. In 
particular, it was shown that the asymmetry can be en- 
hanced by restricting to the small qx region Q. The 
top quark charge asymmetry has received much atten- 
tion recently, due to the deviation from the SM observed 
at the Tevatron 0, @ . Many NP models have been pro- 
posed to explain this discrepancy (see, e.g., Q and ref- 
erences therein). Studying the qr-dependent asymmetry 
will help to clarify which model is the correct one. Simi- 
lar to the asymmetry, it has been shown recently Q that 
vetoing the tt transverse momentum can enhance the sen- 
sitivity of the invariant mass distribution to the effects of 
NP which couples mainly to quarks. This finding makes 
the small qx region even more important. 

Making precise predictions for the small qr region, 
however, is theoretically challenging. As is well-known 
in the case of Drell-Yan and Higgs production, soft and 



collinear gluon emissions give rise to large logarithms of 
the form \n(q T /Q 2 ) at each order in perturbation theory, 
where Q ^> qr is a typical hard scale of the process. The 
fixed-order predictions are therefore not reliable in this 
region. For the case of Drell-Yan and Higgs, the method 
to deal with this problem is the so-called Collins-Soper- 
Sterman (CSS) formalism |{|, in which the large loga- 
rithms can be resummed to all orders in the strong cou- 
pling a s . For tt production, on the other hand, the CSS 
formalism can not be directly applied due to gluon emis- 
sions from the top quarks in the final state. Therefore, for 
observables sensitive to the small qx region in tt produc- 
tion, current experimental groups usually rely on parton 
shower (PS) programs, which only achieves resummation 
at the leading logarithmic (LL) level. Ref. [l(| attempted 
an next-to-leading logarithmic (NLL) resummation by 
modifying the CSS formalism. However, they did not 
consider color mixing between singlet and octet final- 
states, and they missed the contributions from initial- 
final gluon exchange. 

In this Letter, we develop a framework for qx resumma- 
tion in tt production based on the soft-collinear effective 
theory (SCET) [nj]. The framework is built upon the 
works in [l^, EH , which systematically resums the large 
logarithms to arbitrary accuracy. A novel feature of our 
framework is the appearance of a transverse soft func- 
tion matrix, which describes color exchange among the 
initial state and final state particles. Using the available 
ingredients, we perform the resummation at the next-to- 
next-to-leading logarithmic (NNLL) accuracy. 

We consider the process Ni(Pi) + N 2 (P2) -> t(p 3 ) + 
i(p4) + X. We denote the transverse momentum of the tt 
pair as qx- In the small qx region, the differential cross 
section can be written as 
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where s is the collider energy, M and y are the invariant 
mass and the rapidity of the top-quark pair, 9 is the 
scattering angle between p 3 and Pi, fit = \A — 4m 2 /M 2 , 
6,2 = V^ e±a . wi th r = (M 2 + g|)/s. We also define 



Pi = £\P\ , P2 
h = (Pi -P3) 2 



m 2 , 



hi 



s = (pi +p 2 ) 2 , 
(pa -P3) 2 ~ m 2 



The resummed formula for the partonic function Gfi^ ab 
can be written as 



Cij^ ab (zi,Z2,q T ,M,cos0,mt,iJ.) = ^ 



exp [gi(r]i,L±,a s )] 



dbb Jo(6<zt) 



x Tr 
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where r/ z = (C^/tt) ln(Af 2 /^ 2 ) with C q = C F = 4/3 
and C g = C A = 3, L± = ln(&V/& 2 ) with b a = 2e^ E , 
Jo is the zeroth order Bessel function. are the hard 
functions, which are matrices in color space, as indicated 
by the boldface letter. They are the same as in threshold 
resummation and can be found in The functions gi 
and Jj/ are related to the transverse PDFs introduced in 
whose definition and explicit NLO expressions can 
be found in TH, 14 1 . We have recalculated these functions 
at NLO and the results read 
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where T = 4Cj, = — 3CV and 7$ = — (3 = 
— 11Ca/3 + AnfTp/S- Here and below, the perturbative 
functions are expanded in terms of a s /(47r). The split- 
ting functions Pi a {z) are normalized such that Pqq(z) = 
2Cf[(1 + z 2 )/(l — z)] + 1 and the remainder functions 
Ri/ a {z) are given by 
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The functions I' g , a originate from the second Lorentz 
structure of the gluon transverse PDF They start 

at 0(a s ) and do not contribute to the NNLL accuracy. 

The major difference of our resummation formula @ 
from the Drell-Yan or Higgs case is the appearance of the 
transverse soft functions Sg, which are defined as 



S il (L X ,M,co S 9,m u n) = j i Yl J ^d 2 q x e ib ^ 
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where the operator P± acts on the soft final state X s 
giving its transverse momentum, d q i g ^ — 3(8), and Y a 
are soft Wilson lines along the directions of partons a = 
i,i,t,i. The angle <f> t is the azimuthal angle of the top 
quark in the transverse plane. (One may define a soft 
function which is exclusive in (f> t , but the result will be 
much more complicated.) So far we have been working in 
the color space formalism 16[. For actual computations, 
it is more convenient to introduce a color basis, for which 
we adopt the one used in [lj]. In this basis, the LO soft 
functions are given by 
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3. At NLO, the bare soft functions can be 
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where wjj are color matrices, whose explicit expressions 
can be found in [l[ . The integrals Ijk are given by 
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where [dk] = d d k (2tt) 5(k 2 ) 8(k°), and Vj are dimension- 
less vectors along the directions of momenta pj , chosen as 
«i = n= (1,0,0,1), v 2 =n = (1,0,0-1), vj = v% = 1. 
The above integrals contain singularities which are not 
regularized by dimensional regularization. We therefore 
introduce a regularization factor (v/k + ) a following fl7j j . 
where k + = n ■ k, v is an unphysical scale. The results 
for the non-vanishing integrals are then 
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where x s — (1 — /3t)/(l + /3t). The remaining non-zero 
integrals can be obtained by In — Ii3(ti — > u\), I24 — 
123(111 — > t\), I44 = 133. Together with the relations 
among matrices, we find 
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It is worth noting that although the individual integrals 
contain poles in a, these divergences cancel in the final 
soft function, along with the dependence on the unphys- 
ical scale v. The remaining divergences in e are renor- 
malized in the MS scheme. 

Given the renormalization group equations (RGEs) 
satisfied by the hard functions and the transverse PDFs, 
it is straightforward to derive the ones for the soft func- 
tions. We find 



SiM = ~l^(a s )S u (ri - S a (rijMa s ), (10) 
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with 7® = — 2j l l, where jq enter the RGEs of the 
hard functions and can be found in Interestingly, 
the RGEs not only determines the L± dependence of the 
soft functions, but completely fixed them at NLO due to 
the fact that the scale independent terms at this order 
vanish. It would be interesting to see whether this fact 
holds at higher orders in a s . 

Given the resummed formula ([2]), it is important to 
check whether its fixed-order expansion agrees with the 
exact results in the small qx region. To this end we ex- 
pand Eq. fl2J to 0(a s ) and plug it into Eq. (p}. The 
results can be written as 
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where the terms involving the w matrices originate from 
the soft function. 

Eq. (jllj) captures the leading singular terms at order 
a s in the limit qx — >■ 0, which can be compared to the 
exact result in the small qx region. We show in Fig. Q] 
the result from Eq. (fTTl) and the exact result calculated 
using MCFM [l8j . To illustrate the effect of the new soft 
functions, we also show in the plot the result without the 
contributions from the soft functions. As can be seen 
there, only when including the soft function contribu- 
tions, the leading singular terms can reproduce the exact 
result, demonstrating the validity of our formalism. It is 
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FIG. 1. Comparison of the leading singular and the exact 
0(a s ) distributions in the small qr region. Leading singular 
terms with (dashed-dotted line) and without (dashed line) the 
soft function contributions are presented. 



worth pointing out that our Eq. (I12p is in contradition 
with corresponding formulas in [10(. 

To resum the large logarithms, we choose the default 
hard scale as /ih — m t , and evolve the hard function to 
a low scale /i, where the transverse PDFs and the soft 
functions are evaluated. We follow the choice of \i pro- 
posed in [3], [ii = q* + qx for i = q,g, where q* is 
determined by q* = M exp(-2ir/ (r z Q a s (q*))). We also 
adopt the modified power counting such that ol s L\ is 
counted as 0(1). Note that since M > 2m t w 345 GeV, 



q* > 3.28 GeV which is considerably larger than that for 
Z-boson production, where q* « 1.88 GeV. We therefore 
expect much weaker dependence on non-perturbative ef- 
fects in ti production, down to qx = 0. 
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FIG. 2. Resummed predictions for the qx distribution at NLL 
(green band) and NNLL (black band). Also shown are the 
predictions of PDWHEG and MCFM. 

Fig.[5]shows the resummed qx distributions at the NLL 
and NNLL accuracy for ti production at the LHC with 
V 7 * = 7 TeV. Here and below we set m t = 172.5 GeV and 
use MSTW2008NNLO PDFs [H. Uncertainties of the 
theoretical predictions are estimated by varying indepen- 
dently the common scale fj, and the hard scale [ih by a 
factor of two around their central values. It's clear from 
Fig. [5] that the NNLL prediction exhibits significantly 
smaller scale uncertainties, compared with the NLL one. 
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As shown in Fig. [2] the fixed-order prediction from MCFM 
is not reliable when qt is small, while the NLO+PS pre- 
diction of POWHEG [20] is in good agreement with our 
NNLL resummed distribution. It should be noted that 
the POWHEG prediction exhibits a much larger scale de- 
pendence than the NNLL result, which is not shown in 
the plot. 
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FIG. 3. Comparison of NNLL resummed prediction (blue 
band) for the normalized qr distribution with the experimen- 
tal data from the CMS collaboration. 
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FIG. 4. The top quark charge asymmetry as a function of qr- 
The Pythia and MC0NL0 curves are extracted from [4]. 

In Fig. [3] we show our NNLL resummed prediction for 
the normalized qr distribution, together with the ex- 
perimental data from the CMS collaboration Q, using 
an integrated luminosity of 1.14 fb _1 at the LHC with 
yG = 7 TeV. In this plot a non-perturbative factor is 
included for the gg-channel but not for the gg-channel, 
where the non-perturbative effects are believed to be neg- 
ligible. The experimental data shows good agreement 
with our resummed prediction. 

We finally turn to the gT-dependent top quark charge 
asymmetry Afb ■ This quantity is of substantial interest 
because it will provide new hints for the puzzle of large 
deviation in A?b observed at the Tevatron. In QCD, the 
asymmetry starts at NLO, however, it was found that an 
LO parton shower program like Pythia can exhibit non- 
zero Afb ■ As was explained in 2l| , this is due to the fact 
that in the hard process qq — ¥ ti, color coherence of par- 



momentum when the top goes backwards. In our resum- 
mation formalism, this color coherence is accounted for 
by the soft function S q q, whose dependence on t\ and u\ 
is asymmetric. In Fig. 01 we present our resummed pre- 
diction for this observable, together with predictions from 
MC0NL0 and Pythia extracted from Interestingly, our 
NNLL resummed prediction shows very good agreement 
with the NLO+PS program MC0NL0. In particular, they 
predict the same cross-over at qr ~ 25 GeV. 

In conclusion, for the first time, we have presented 
a resummation framework for the transverse-momentum 
spectrum of top-quark pairs at hadron collider, valid up 
to arbitrary logarithmic accuracy. Compared with Drcll- 
Yan and Higgs production, a new ingredient in our for- 
malism is the introduction of the transverse soft function 
matrices, which describe the soft gluon effects associated 
with final-state radiations. We have explicitly shown that 
when expanded to 0(a s ), our resummation formula re- 
produces precisely the fixed-order prediction from MCFM 
at small qx- We have carried out the resummation at 
NNLL accuracy. Our results agree quite well with those 
from parton shower programs and with the CMS mea- 
surement, while exhibiting a small scale dependence. We 
have also examined the gr-dependent top quark charge 
asymmetry, which could help clarifying the large devia- 
tion from the SM observed at the Tevatron. Our formal- 
ism can also be applied to the bb, cc production, as well 
as the production of colored supersymmetric partners. 
With the NNLO soft function which may be calculated in 
the future, our work provides a new subtraction method 
for computing the tt differential cross sections at NNLO, 
following the qx subtraction method of [22j. Finally, it 
is interesting to incorporate the decays of the top quark 
into our framework in a way similar to [23], which we 
leave for future works. 
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